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ABSTRACT 

Nearby spiral galaxies show an extremely tight correlation between tracers of molecular hydrogen 
(H2) in the interstellar medium (ISM) and tracers of recent star formation, but it is unclear whether 
this correlation is fundamental or accidental. In the galaxies that have been surveyed to date, H2 
resides predominantly in gravitationally bound clouds cooled by carbon monoxide (CO) molecules, 
but in galaxies of low metal content the correlations between bound clouds, CO, and H2 break down, 
and it is unclear if the star formation rate will then correlate with H2 or with some other quantity. 
Here we show that star formation will continue to follow H2 independent of mctallicity. This is not 
because H2 is directly important for cooling, but instead because the transition from predominantly 
atomic hydrogen (H i) to H2 occurs under the same conditions as a dramatic drop in gas temperature 
and Bonnor-Ebert mass that destabilizes clouds and initiates collapse. We use this model to compute 
how star formation rate will correlate with total gas mass, with mass of gas where the hydrogen is 
H2, and with mass of gas where the carbon is CO in galaxies of varying metallicity, and show that 
preliminary observations match the trend we predict. 

Subject headings: galaxies: star formation — ISM: atoms — ISM: clouds — ISM: molecules — stars: 
formation 



1. INTRODUCTION 

Spiral galaxies near the Milky Way show strong cor- 
relation between the surface density of star formation 
and the surface density of molecular gas, but only a very 
weak correlation betw een star formation and the surface 
density of a tomic gas f Kennicutt et al.ir2007t iBigiel et"al] 
l2008HLero v ct al. 2008). However, the origin of this be- 
havior is highly deba ted. In many models for the star 
formation rate OTanll2000l: ILi et al.fe ool ISilk fc NormanI 
120091: IDobbs fc Pringlell2009f) . the chemical state of the 
gas is assumed to be dynamically unimportant. In these 
models stars form wherever galactic-scale gravitational 
instabilities dictate that gas collects into gravitationally 
bound structures, without regard to whether those struc- 
tures consist of atomic or molecular gas. The observed 
correlation between molecules and star formation occurs 
only because molecules form preferentially in the same 
bound structures where stars do. 

In other mo dels (jRobertson fc Kravtsovl I200S : 
Gnedin et all 120091: 'Krumh olz et all l2009di 



Pelupessv fc Pao adopouloa i2009. : iGnedin fc Kravtsovl 
2OIOI : iFu et al. 20i3), however, the underlying assump- 



tion is that star formation follows the chemical transition 
between atomic hydrogen (H i) and molecular hydrogen 
(H2) as the dominant gas phase, or the transition from 
ionized carbon (C 11) to carbon monoxide (CO) as the 
primary coolant, independent of the global dynamics. 
Closely related to these are models in which the locations 
where stars can form is determined by the gas tempera- 
ture, or by the interplay betw een gas temperatu r e and 
gravitational stability ((Elmegreen fc Parravanol 119941 : 
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lSchavdl2004lOstriker et~alll2010t ). It is not entirely clear 
how this is corre l ated w ith the chemical composition, 
although 'Schay^ (|2004[ ) obtain the intriguing result, 
which anticipates to some extent the one we report here, 
that the formation of a cold phase is closely associated 
with the H I to H2 transition. However, he did not 
consider how this relates to CO, the molecule usually 
observed in place of H2 . 

All of these models produce similar results for nearby 
spiral galaxies, where the majority of gravitationally- 
bound structures are also molecular and cooled by CO. 
However, whether star formation follows global dynam- 
ics, the H I to H2 transition, or the C 11 to CO tran- 
sition makes a very large difference for how star for- 
mation behaves in regions with substantially sub-Solar 
metal content, since the chemical states of the hydrogen 
and the c arbon are sensitive to gas me t allicity in differ- 
ent wa ys (Ivan Dishoeck fc Black 11988 1: iKrumhol z et al.l 
I2009d : iGlover fc Mac Lowl 12010 : IWblfirc et al. 2010H. 
while the prevalence of galactic-scale instabilities is very 
insensitive to metallicity. Moreover, any successful model 
must also be able to relate the star formation rate to ob- 
servable quantities, which include H i and CO, but not 
(generally) the total mass of H2. 

At present observations possess limited power to dis- 
criminate between the models, although the limited 
data available for nearby low-metallicity dwarf galax- 
ies such as the Small Magellanic Cloud hold out the 
promise of being able to distinguish whether metallicity 
affects star formation. However, the differences between 
the models have profound implications for our under- 
standing of star formation in the high redshift universe, 
where average metallicities are much lower. For exam- 
ple, the difficulty of forming molecules in metal-poor gas 
has been proposed a s an explanation for the observed 
(jWolfe fc ChenI 120061 : IWild et al.ll2007l ) lack of star for- 
mation in damped Lym an a systems (jKrumholz et all 
,2009ai : IGnedin fc Kravtsov 2010) . but this explanation 
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is only viable if transitions to molecular gas are neces- 
sary for star formation. 

Here we investigate whether the chemical state of the 
gas is relevant for star formation by studying how the 
Bonnor-Ebert mass and the chemical state vary in spheri- 
cal interstellar clouds of varying volume and column den- 
sities. We show that, for such clouds, the H I to H2 
transition is strongly associated with a dramatic drop 
in the temperature and Bonnor-Ebert mass within the 
cloud, across a very wide range of metallicities and en- 
vironments. In Section [2] we describe our computation 
method, in Section |3] we report our results, and in Sec- 
tion 21 we discuss the implications of our results for ob- 
servations. Finally, we summarize in Section [5] 

2. COMPUTATION METHOD 

To minimize clutter we summarize all the physical pa- 
rameters that enter our model in Table (TJ and we dis- 
cuss these choices in Appendix |Al We also show in that 
Appendix that our results are robust against changes in 
these parameters. In this section we limit ourselves to a 
discussion of our calculation method. 

2.1. Physical Model 

Consider a spherical cloud of mean volume density nn 
H nuclei cm"'^ and mean column density Nn H nuclei 
cm~^, mixed with dust with an absorption cross section 
per H nucleus aj. to photons near 1000 A. The mean 
absorption optical depth to UV photons is r = N^ad, 
and the corresponding mean visual extinction at optical 
wavelengths is Ay- The outer parts of the cloud are ex- 
posed to the ultraviolet interstellar radiation field (ISRF) 
of the galaxy, and this keeps the hydrogen and carbon in 
the outer parts of the cloud predominantly in the form of 
H I and C 11. If the column density is sufficiently large, 
dust and the small population of hydrogen molecules in 
the predominantly atomic layer will absorb the ISRF, al- 
lowing a transition from H i to H2 and C 11 to CO as 
the main repositories of hydrogen and carbon. We as- 
sume that the gas is neutral, and that there are very few 
ionizing photons present. 

In a real interstellar cloud, even if the density were uni- 
form, the temperature and chemical composition would 
not be. Surface layers exposed to direct UV radiation 
would be warmer than the shielded interior, and the hy- 
drogen there would be primarily H i rather than H2. 
Thus there is no single cloud temperature or chemical 
composition we can compute, and determining the full 
temperature and chemical profile even when the density 
profile is given in advance requires solving a photodisso- 
ciation region model in which one both simultaneously 
determines temperature and the chemical state, includ- 
ing line radiative transfer throughout the cloud. While 
there are many examples of such models in the literature, 
they are too computationally expensive to allow a wide 
search of parameter space, and they still rely on arbi- 
trarily chosen density distributions. Instead, our goal is 
to estimate a single characteristic temperature and for a 
cloud, coupled to a simple description of its chemistry. 
To this end we will describe clouds' chemical composition 
with single numbers giving the mass fraction in a given 
chemical phase, and we will approximate the cloud inte- 
rior as a uniform region of volume density rifj and mean 
column density A'h, consisting of gas at temperature Tg 



uniformly mixed with dust of temperature T^. For the 
purposes of computing the temperature, we will take the 
chemical composition to be uniform, with the hydrogen 
all as H I or H2, and the carbon all as C 11 or CO, despite 
the fact that every II2 cloud has an outer H i shell, and 
every CO cloud has an outer C 11 shell. 

2.2. Chemical State 

To compute the H i to H2 transi tion in our mod el cloud, 
we r ely on the models of Krum holz et all ()2008l [2009c), 
and iMcKee k Krumholz (^2010) (collectively the KMT 
model hereafter). In these models the fraction of the 
hydrogen mass in the H2-dominated region is given by 



/h2 



1 - 



1 + 0.25s 



where 



ln(l + 0.6x + 0.01x^) 



0.6r 



X = 71 



G'n 



nil /cm ^ 



(1) 

(2) 
(3) 



and Gq is the strength of the ISRF normalized to its 
value in the Solar vicinity (which corresponds to a free- 
space H2 dissociation rate of 5.4 x 10~^^ s~^). Note 
that, unlike some of the other approximate expres- 
sions derived by KMT, this one applies independent of 
whether the gas is warm or cold neutral medium. It 
should also be noted that this model assumes chemi- 
cal equilibrium, which is not necessarily the case for H2 
tPelupessv & Papa dopoulos 2009: „M ac Low fc Glovei) 
120101 )^ However, Krumholz fc GnedinI |20lH ) perform 



a detailed comparison between the equilibrium approx- 
imation and a set of galaxy simulations including full 
time-dependent chemistry and radiation. They find that 
the approximation is very accurate at metallicities of 
Z' > 0.01, where Z' is the metallicity normalized to 
the Milky Way value, and we can therefore use it safely. 
We note that equation ([T]) also agree s extremely well 
with o bservations both of local galaxies ([Krumholz et al.l 
l2009cD and of damped Lyma n a systems at high red- 
shift ([Krumholz et al.l l2009b[) . We therefore conclude 
that, on the galactic scales relevant to this work (as 
opposed to the isolated periodic boxes considered by 
iMac Low Sz Gloverll2010l ) this approximation is valid. 

The C II to CO transition is predominantly governed 
by dust extinction, and thus is sensitive primarily to the 
cloud optical depth. At high optical depths wherever the 
hydrogen is H2 the carbon is also CO, while at low optical 
depths there can be significant regions of H2 where the 
carbon is primar ily C il, because the H2 self-shields while 
the CO cannot ([Glover fc Mac Lowl[20Tot iWolfire et all 
[2OT0f l. Recent semi-analytic work indicates that the ratio 
of mass where the carbon is CO to total cloud mass is 
(IWolfire et al.l[201Q) 



-4.0 0.53-0.045 In " _j -0.097 In Z' /Av 

fcO = /HaC V "H/cm J 

... (4) 

Alternately, an em pirical fit to numerical sim ulations for 
the same quantity (I Glover fc Mac Lowl[20T0l) gives very 
similar results. 
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TABLE 1 

Fiducial parameter choices 



Parameter 



Value 



Meaning 



Reference 



0.01(T/10 K)^Z' cm2 g" 



"gd {H2) 


3.2 X 10-34^' erg cm^ 




Qgd (H I) 


1.0 X IQ-^^Z' erg cm^ 




9CR (H2) 


12.25 eV 




"JCR (H l) 


6.5 eV 




c 


2 X 10-1'^ Z' s-i 




Trad 


8 K 




G'g 


1 






1.6 X 10-* Z' 






1 X 10~2iz' cm" 


2 


Av/Nii 


4.0 X 10-22^' mag 


cm^ 


OPR 


0.25 






2 km s-l 




Note. — References to AppendixlAI mean 

2.3. Thermal State 



^ IR dust opacity 

Dust-gas heat exchange rate 
^^'^ Dust-gas heat exchange rate 
Energy / CR ionization 
Energy / CR ionization 
CR ionization rate 
Ambient radiation temperature 
UV radiation intensity 
C II / CO abundance 
UV dust opacity per H nucleus 
Visual extinction per H column 
Ratio of ortho- to para-H2 
Cloud velocity dispersion 



Lesaffre et al. (2005) 
Goldsmith (20011 
Appendix 
Appendix lAI 
Dalearno fc McCrav (1972) 
^^^^^^ppendix^y 
Appendix ^ 
Appendix \K\ 
Sofia et al. (2004) 
AppendixlAI 
Appendix [XI 
Neufeld et a l. (2006) 
Krumholz fc Th ompson (2007) 



The temperatures of gas and dust are set by the 
condition of thermal equihbriu m, followin g a m ethod 
that combine s elem ents from iGoldsmithI (|2001[ ) and 
iLesaffre et al.l (|2005[ ). We consider heating of the gas 
by the grain photoelectric effect at a rate per H nucleus 
Fpe, heating by cosmic rays at a rate Fcr, and cooling 
via atomic and molecular line emission at a rate Aiinc- 
Dust is heated via an external radiation field at a rate 
Fdust, and cools via thermal emission at a rate Adust- 
Finally, energy flows from the dust to the gas at a rate 
VPgd, which may be negative if the gas is hotter than 
the dust. The condition for thermal balance is therefore 
that the temperatures Tg and Td simultaneously satisfy 
the equations 



Tdust ~ Adust 



*gd = 



(5) 
(6) 



Note that we do not include photoionization heating. 
While this i s important at low column densities (e.g. 
iSchavd [20041 ) ■ the lowest column density clouds we will 
consider in this work have optical depths of ^ 10^ to ion- 
izing photons, and thus photoionization heating is unim- 
portant for them except very close to their surfaces. 

To solve this equation we must compute the tempera- 
ture dependence of each heating and cooling term. We 
assume that the clouds in question are optically thin to 
far-infrared radiation, so we need not consider trapping 
of the dust radiation field. Thus the dust cooling rate is 



A, 



dust 



(7) 



where c is the speed of light, a is the radiation constant, 
n{Td) is the temperature-dependent dust specific opacity 
and /iH is the mean mass per H nucleus. Dust heating is 
more complex, as mentioned above, since the external in- 
terstellar ultraviolet field dominates at cloud edges and 
the local re-radiated infrared field dominates at cloud 
centers. Since we are mostly interested in what happens 
when clouds are shielded by at least a few magnitudes of 
extinction in the ultraviolet, we choose to adopt the IR- 
dominated case. In practice this choice makes very little 
difference, since the dust gas coupling is a very minor 
effect at the densities where we focus our attention. We 
characterize the IR radiation field by an effective tem- 



perature Trad, so 

Tdust 



'«(T'rad)/^HCaT4d- 



(8) 



Finally, for the dust-gas energy exchange term w e adopt 
the approximate exchange rate (|Goldsmithll200lD 



gd 



where ofgd is a coupling constant. 
For the gas, the cosmic ray heating rate is 

Tcr = C^CR s"^ 



(9) 



(10) 



where C is the cosmic ray primary ionization rate and 
qcK is the thermal energy increase per primary cosmic 
ray ionization. For grain photoelectric heating, we must 
choose a suitable mean heating rate, since extinction 
through the cloud causes the heating rate to drop sub- 
stantially as we move toward the interior. For simplicity, 
and since we are concerned more with cloud interiors 
than surfaces, we consider the heating rate to be attenu- 
ated by half the mean extinction of the cloud. With this 
approximation, we have 



Fpe = 4.0 X 10" 



-26 



(11) 



where ad is the dust cross section per H nucleus to the 
UV photons that dominate grain photoelectric heating. 

The final term in the thermal balance equations ([5]) 
and ([B]) is Aiiuc, the line cooling rate. In each case we 
consider only a single coolant species: C 11 or CO. This 
is a reasonable approximation because, at the low tem- 
peratures where star formation occurs, one of these two 
will dominate. To compute the cooling rate we must 
compute the populations of th e various levels, and we 
do so following the method of jKrumholz fc ThompsonI 
(j2007[) . Let fi be the fraction of a given species, C li or 
CO, in the ith state, and let Ei be the energy of that 
state. Transitions between states i and j occur due to 
spontaneous emission, at a rate given by the Einstein 
coefficient Aij (which is zero if Ei < Ej), and due to 
collisions with with rate coefficient 7^.3, where S is the 
species of the collision partner. Note that the collision 
rate coefficients depend on the gas temperature. We do 
not consider stimulated emission or absorption due to an 
external radiation field, since neither is significant for the 
cooling lines of C 11 or CO. 
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In the escape probability formalism we approximate 
that each spontaneously emitted photon produced by an 
atom transiting from level i to level j has a probability 
of escaping; photons that do not escape are re-absorbed 
on the spot within the cloud, yielding no net change in 
the level populations. With these definitions and approx- 
imations, in statistical equilibrium the level populations 
are given implicitly by 

fi X! (l^V^ij + nmji],m + ?T-He7u,He 

3 

+ «p-H27y,p-H2 + 'T-o-H27ij\o-H2) = 



^p-H27iJ,p-H2 + '^o-H27'iJ,o-H2) fi 



(12) 



where the sums run over all quantum states, and nni, 
np_H2, "■0-H2J and 'T'He are the number densities of H i, 
para-H2, ortho-H2, and He, respectively. The left hand 
side of this equation represents the rate of transitions 
out of state i to all other states j, while the right hand 
side represents the rate of transitions into state i from 
all other states j. The escape probabilities themselves 
depend on the level population. If we let X be the abun- 
dance of the species in question relative to H nuclei, and 
a be the one-dimensional gas velocity dispersion, then 
we have ([Krumholz fc ThomDSon|[2007[ ) 



where 



A, 



3 j 



1 



1 + Q.bT,, 



gj 16(27r)3/2a 



(13) 



(14) 



Qi and gj are the statistical weights of levels i and j, and 
Xij is the wavelength of a photon associated with the 
transition from level i to level j. 

Equations - ^1)1 constitute a complete system of 
algebraic equations for the level populations fi. Given 
the solution to these equations, the line cooling rate per 
H nucleus is then given by 



Aline — PijAijfihvij/nji, 



(15) 



where Vij is the frequency of a photon emitted in a tran- 
sition from state i to state j. 

We have now written down a complete set of equa- 
tions to determine the gas and radiation temperatures. 
To solve them we use a double-iteration method. We 
first select trial values of Td and Tg, and then we com- 
pute all the heating and cooling rates F, A, and ^. In 
order to compute A^ne , we must solve for the equilibrium 
level populations fi by solving equations - ([T^ . We 
do so by fixing Tg (and thus all the collision rate coef- 
ficien ts) and applying a Broy den's method (Pres s et al.l 
Il992| ). Once we have determined Aunc, we check if the 
thermal balance equations ([5]) and © are satisfied to 
within a specified tolerance. If not, we iteratively update 
Td and Tg using Newton's method to update Td and Tg 
until equations ([5]) and ([6]) are satisfied to the required 
tolerance. 




Fig. 1. — Color indicates gas temperature as a function of density 
njj and visual extinction Ay for metallicities of Z' = 0.05, 0.2, 
and 1.0 (left to right columns; prime indicates relative to the Solar 
neighborhood value) and cooling by C II plus H I (top row), C II plus 
H2 (middle row), and CO plus H2 (bottom row) using our fiducial 
parameters Gg = 1, C-17 = C/10~^^ = 2Z' (upper left corner). 
Color contours run from log(r/K) = 0.5 to log(T/K) = 2.5 in steps 
of 0.1, with white representing log(T/K) > 2.5. The black solid 
lines are contours of fn^ = 0.1, 0.5, and 0.9, from left to right. 
The gray dashed lines are contours of fco = 0-1 E'nd 0.5 from 
left to right; the /co = 0.9 contour lies off the plot to the right. 
Finally, the hatching pattern marks regions where the temperature 
and chemistry calculations are not fully consistent (e.g. because 
the chemistry calculation indicates that a cloud of that Ay and 
n should be dominated by H2, but the temperature calculation 
assumes the hydrogen is H l). 



3. RESULTS 

We use the procedure outhned in Section [2] to com- 
pute the chemical state and temperature for a grid of 
clouds of varying nn and Ay ■ Figure [1] shows the equi- 
librium temperature as a function of ripj and Ay , overlaid 
with contours of H2 and CO fraction, computed for the 
three possible chemical compositions and gas metallici- 
ties of Z' — 0.05, 0.2, and 1.0. Recall that the chemical 
and thermal computations are decoupled, so the temper- 
atures for C II plus H I cooling should only be considered 
rchable to the left of fco — 0.5 and /h2 = 0.5 contours, 
indicating the C 11 to CO and H i to H2 transitions, re- 
spectively. Similarly, the CO plus II2 temperatures are 
only reliable to the right of both of these curves, while 
the C II plus II2 temperatures are reliable in the region 
between them. 

Regardless of these limits on the regions of applica- 
bility, the striking result from these plots is that, inde- 
pendent of metallicity or chemical composition, there is 
a dramatic drop in temperature from hundreds of K to 
^ 10 K as one moves from lower left (low density, low 
Ay) to upper right (high density, high Ay), and that 
the contours of constant temperature align remarkably 
well with contours of constant /hj • This result is not 
surprising, because temperature and II2 fraction depend 
on density and extinction in very similar ways. Both the 
heating and II2 dissociation rate are proportional to the 
the exponential of minus the visual extinction, and both 
the cooling and H2 formation rates are proportional to 
the square of the volume density. In contrast, contours of 
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Fig. 2. — Same as Figure [T] but with color indicating Bonnor- 
Ebert mass instead of temperature. Color contours run from 
log(MBE/M0) = 0.5 to log(MBE/A^0) = 3.5 in steps of 0.25, with 
white representing log(MBE/A^Q) > 3.5. 




3.0 0.2 0.4 0.6 O.S 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0. 



Fig. 3. — Contours are the same as in Fig. 2, but the x axis 
shows H2 fraction rather than Ay. Note that not all of the grid 
is colored, because our grid in nu and Ay does not cover the full 
range in njjj and /jjg shown. 

constant CO align much less well with temperature con- 
tours. The CO fraction depends primarily on Ay, with 
only a very weak density dependence. Furthermore, at 
Av high enough for the CO fraction to reach 50%, the 
grain photoelectric effect has been shut off so thoroughly 
that the temperature is insensitive to further increases 
in Ay. 

Figures [2] and [3] show how the Bonnor-Ebert mass, the 
largest mass that can be supported against collapse by 
thermal pressure, changes with Ay and with H2 fraction 
as a result of this temperature change. This mass is 



BE 



1.18 



= 1.18 



1/2 



(16) 



where Cg is the isothermal sound speed, T is the temper- 
ature, and /LtH and n are the mean mass per H nucleus 



and the mean particle mass, respectively. For Milky Way 
helium abundance, the former is = 2.3 x 10"^''' g 
regardless of chemical composition, while the latter is 
2.1 X 10-24 g for H I and 3.8 x lO'^"* g for H2. As the 
plot shows, there is a radical drop in the Bonnor-Ebert 
mass from thousands of Mq to a few M© as one moves 
from low density and extinction to high density and ex- 
tinction, and, as with the temperature, there is a very 
strong correlation between Bonnor-Ebert mass and H2 
fraction, while there is relatively little correlation with 
CO fraction. Figure 3 clearly shows that small Bonnor- 
Ebert masses are found exclusively in clouds with high 
H2 fractions. 

Our results imply that, in gas whose conditions are 
such that the hydrogen is mostly H i, structures with 
masses of ~ 1000 M© or less will be stabilized against 
collapse by thermal pressure. As a result, we conclude 
that star formation in such environments is extremely un- 
likely, in part because turbulence is extremely unlikely to 
generate fragments of such high ma sses that could then 
collapse (jPadoan fc Nordlundl 120021 ) . Conversely, in gas 
where the hydrogen is mostly H2, the mass that can be 
stabilized by thermal pressure is 2— 3 orders of magnitude 
smaller, turbulence generates numerous fragments capa- 
ble of collapsing, and star formation is far more likely. It 
is important to note that the drop in Bonnor-Ebert mass 
and loss of stability is not caused by the H i to H2 tran- 
sition, it is simply very well correlated with it because 
the temperature and the H2 fraction are determined by 
very similar combinations of extinction and density. We 
therefore conclude that star formation should correlate 
well with H2, simply because H2 is a good tracer of re- 
gions where thermal pressures are low enough to allow 
star formation. 

4. OBSERVATIONAL CONSEQUENCES 
4.1. Predictions for Observable Quantities 

We can use our result that star formation correlates 
with H2 to predict an approximate correlation between 
star formation rates and the masses of all gas, H2 gas, 
and CO g function of galaxy metallicity and sur- 

face density. Consider a portion of a galaxy with a mean 
surface density (averaged over the ~kpc scales acces- 
sible to current olsservations). Within this region a frac- 
tion of the gas, given by equation ([T]), is in the form 
of star-forming H2 clouds, and within these a smaller 
fraction fco of the gas, given by equation Q, has most 
of its carbon in CO molecules. 

We evaluate fa, follo wing the methods outlined in 
iKrumholz et al.l (|2009cl ldl). First, we must adopt a 
clumping factor c to scale from the surface densities 
of individual atomic molecular complexes on ~ 100 pc 
scales to the ^^kpc scales accessible to current obser- 
vations. This is necessary because the mean surface 
density averaged over a kpc-scale region that we ob- 
serve is lower than the surface densities of the ^ 100 
pc-sized giant molecular clouds, but it is the latter 
rather than the former that determines the atomic to 
molecular ratio in these clouds. Based on a combina- 
tion of theoreti c al arg uments and fits to observation, 
IKrumholz et al.l ()2009c[ ) adopt c = 5, and we do so here, 
so T = 5Eg(Trf//XH- Second, we adopt a characteristic 
value of G'o/n = 0.044(1 + 5.lZ'°-^^^)/AA expected for 
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a two-phase atomic medium ()Krumholz at al.ll2009cl ldl). 
(For a given value of Gn, this corresponds to assuming 
that clouds on Figures [J - [3] form a horizontal sequence 
at a fixed n value that depends on G'q and Z'.) To eval- 
uate /co s-nd the star formation rate, we note that star- 
forming H2 clouds appear to develop column densities 
Nh 7.5 X 10^1 cm"^ 85 M(d pc~^) independent 
of ga, l actic environment (iHeyer et al]|2009l : iBolatto et alj 
120081: iKrumholz eFall 120061). Since the CO-dominated 
regions are the inner parts of these clouds, this implies 
that the column density and Ay that enter the calcula- 
tion of the CO fraction should not be the mean galactic 
surface density or extinction, but instead the fixed Ay 
of star-forming H2 clouds. We adopt this value for Ay 
in equation (|4]). 
The corresponding star formation rate is ex- 

ected on the oretical grounds to be approximately 

Krumholz et al. 2009d) 



2.6 Gyr 

(S3/85M0pc-^) 
(I]g/85M0pc " 



-0.33 



Eg < 85 Mq pc" 



2^0-33 Eg > 85M0PC-2 (l'^) 



In reality, based on the argument we have just made, 
we should compute the star formation rate not based 
directly on the H2 fraction but instead based on the mass 
of cold gas. However, we have already seen that the 
H2 and cold gas fractions are very similar, and equation 
([1} provides a convenient analytic approximation. We 
therefore use it to estimate both and the star-forming 
gas fraction. We have therefore computed, for a portion 
of a galaxy of total gas surface density Eg and metallicity 
Z', the expected surface densities of H2, CO, and star 
formation. 

Figure H] shows our predicted correlation between spe- 
cific star formation rate with respect to the mass of each 

gas constituent, SSFR(co32, total) = E,/E(co,H2,totai), 
and the surface density of that component at a range of 
metallicities. We see that at high surface densities and 
high metallicities the specific star formation rates for to- 
tal gas, H2, and CO are essentially identical, consistent 
with observations. At lower surface densities or metallic- 
ities, though, the star formation rate per unit total gas 
falls. Conversely, at lower metallicity the specific star for- 
mation with respect to CO rises, reflecting the fact that, 
at lower metallicity, the mass of gas where the chemical 
makeup is H2 plus C 11 rises as a fraction of the total H2 
mass. Since star form ation follows H2 , SSFRpo rises as a 
result. We note that iPelupessv &: Papadopoulo^ (|2009l ) 
reached a similar conclusion about SSFRco based on 
simulations that used a star formation recipe that as- 
sumed roughly constant SSFRna- 

4.2. Comparison to Observations 

The most striking feature of Figure U is that it pre- 
dicts very strong variation in specific star formation rate 
with metallicity for CO and total gas, but not for H2. In 
the low metallicity galaxies where we would like to test 
this prediction it is difficult to measure either the total 
gas mass or the H2 mass, because we lack a reliable H2 
tracer except in a few limited cases (see Section 14.31 for 
further discussion). We can, however, measure CO lu- 
minosities and thereby estimate masses of gas that are 
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Fig. 4. — Predicted gaseous specific star formation rate 



SSFR, 



(CO, H2 .total) 



(CO, H2, total) 



CO,H2, total) 



a function of gas surface density at fixed metaUicity. The red fines 
indicate the specific star formation rate considering the total gas 
mass, green lines indicate the specific star formation rate for H2, 
while blue lines indicate the specific star formation rate for gas 
where the carbon is CO rather than C 11. Dashed lines correspond 
to a metallicity Z' = 0.05, dotted lines to Z' = 0.2, and solid lines 
to Z' = 1. Note that the dashed blue line (SSFRco for Z' = 0.05) 
is above the range shown in the plot. 



CO-dominated. As shown in Figure 21 we predict that 
the specific star formation rate for CO should be much 
higher in the lower metallicity systems. 

To test this prediction, we compile galaxy-integrated 
CO luminosities, Lqq, and star formation rates, SFRs, 
from the literature. These measurements cover a wide 
range of metallicity from highly sub-Solar to super-Solar, 
and include both spirals and dwarfs. The data come from 
a variety of sources, summarized in Table [2j Whenever 
they are available, we give preference to Leo derived 
from integrating complete single-dish maps of a galaxy, 
but many of the data still come from sparse sampling. 
For SFRs, we draw heavily f rom the recent s et of galaxy- 
integrated measurements by iCalzetti et al.l ()2010D ; these 
have the advantage of being computed in a uniform way 
and include an IR contribution to correct for extinction. 
Both Lqo and the SFR are luminosity-like quantities, 
so the ratio of the two is independent of distance. We 
also draw metallicities for each target from the literature. 
When avail able, we give pre f erence to the recent com- 
pilation by iMo ustakas et al.l ()2010[) , using the average 
of their two characteri stic m etallicities. For targets not 
studied by Mousta kas et al. . we use values from the com- 



pila tions bv Marble et all (I2010D . ICalzetti et all (|2010t ). 
and lEngelbracht et al.l (|2008[ ) and a variety of literature 
sources. 

To convert observed CO (1 — 0) luminosities to masses 
of gas traced by CO, we adopt a conversion factor 2 x lO^*^ 
Ha molecules cm-V(K km s~ ])-^ (lAbdo et all l2010bl : 
iBlitz et al.l[2007t iDraine et al][2007L iHever et al.ll2009( ). 



so the mass is 



Mco = 4 X 10^ 



^co 



K km s cm^ 



(18) 



where /xh is the mean mass per H nucleus. This is equiv- 
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TABLE 2 

Data compilation 



Galaxy log Leo [K km s— 1 pc^] CO reference logAf, [-^0 yr ^] SFR reference 12 + log{0/H) Metallicity reference 



SMC 


5.2 ± 0.2 


M06 


-1.3 ±0.2 


W04 


8.0 ±0.2 


D84; MAIO 


LMC 


6.5 ±0.1 


F08 


-0.7 ±0.2 


H09 


8.3 ±0.3 


D84; MAIO 


ICIO 


6.3 ±0.2 


L06 


-1.0 ±0.5 


L06 


8.2 ±0.2 


L79; L03 


M33 


7.6 ±0.1 


H04 


0.0 ±0.5 


H04 


8.3 ±0.2 


R08 


IIZw40 


6.2 ± 0.3 


T98 


-0.2 ±0.3 


CIO 


8.1 ±0.3 


E08; CIO 


NGC1569 


5.1 ± 0.3 


T98 


-0.6 ±0.3 


CIO 


8.1 ±0.2 


M97 


NGC2537 


5.5 ± 0.3 


T98 


-1.1 ±0.3 


CIO 


8.4 ±0.3 


MAIO 


NGC4449 


6.8 ± 0.3 


Y95 


-0.5 ±0.3 


CIO 


8.3 ±0.2 


M97 


NGC5253 


5.8 ± 0.3 


T98 


-0.3 ±0.3 


CIO 


8.2 ±0.3 


MAIO 


NGC6822 


5.1 ± 0.2 


197 


-2.0 ±0.3 


CIO 


8.4 ±0.3 


MOlO 


NGC0628 


8.3 ± 0.1 


L09 


-0.2 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC0925 


7.6 ± 0.1 


L09 


-0.3 ±0.2 


CIO 


8.6 ±0.3 


MOlO 


NGC1482 


8.8 ± 0.3 


Y95 


0.5 ±0.2 


CIO 


8.6 ±0.4 


MOlO 


NGC2146 


9.3 ± 0.3 


Y95 


0.9 ±0.2 


CIO 


8.7 ±0.3 


E08; CIO 


NGC2403 


7.1 ± 0.3 


Y95 


-0.4 ±0.2 


CIO 


8.6 ±0.2 


MOlO 


NGC2841 


8.4 ±0.1 


L09 


0.1 ±0.2 


CIO 


9.0 ±0.3 


MOlO 


NGC2782 


9.0 ±0.3 


Y95 


0.7 ±0.2 


CIO 


8.6 ±0.4 


E08; CIO 


NGC2798 


8.8 ±0.3 


Y95 


0.5 ±0.2 


CIO 


8.7 ±0.3 


MOlO 


NGC2976 


7.0 ±0.1 


L09 


-1.0 ±0.2 


CIO 


8.7 ±0.3 


MOlO 


NGC2903 


8.8 ±0.1 


H03 


0.3 ±0.2 


CIO 


8.9 ±0.3 


MAIO 


NGC3034 


8.8 ±0.3 


Y95 


0.9 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC3077 


6.0 ±0.3 


T98 


-1.0 ±0.2 


CIO 


8.6 ±0.3 


MAIO 


NGC3079 


9.4 ±0.3 


Y95 


0.5 ±0.2 


CIO 


8.6 ±0.4 


E08; CIO 


NGC3184 


8.4 ±0.1 


L09 


-0.5 ±0.2 


CIO 


8.9 ±0.3 


MOlO 


NGC3198 


7.9 ± 0.1 


L09 


-0.0 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC3310 


8.2 ± 0.3 


Y95 


0.9 ±0.2 


CIO 


8.2 ±0.4 


E08; CIO 


NGC3351 


8.2 ± 0.1 


L09 


-0.2 ±0.2 


CIO 


8.9 ±0.3 


MOlO 


NGC3368 


8.3 ±0.3 


Y95 


-0.4 ±0.2 


CIO 


9.0 ±0.3 


MAIO 


NGC3521 


8.8 ±0.1 


L09 


0.1 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC3628 


9.2 ±0.3 


Y95 


0.3 ±0.2 


CIO 


9.0 ±0.3 


MAIO 


NGC3627 


9.0 ±0.1 


H03 


0.2 ±0.2 


CIO 


8.7 ±0.3 


MOlO 


NGC3938 


8.5 ±0.1 


H03 


-0.1 ±0.2 


CIO 


8.7 ±0.3 


E08; CIO 


NGC4194 


8.9 ± 0.3 


Y95 


1.1 ±0.2 


CIO 


8.7 ±0.4 


E08; CIO 


NGC4214 


6.1 ± 0.1 


L09 


-1.0 ±0.2 


CIO 


8.2 ±0.2 


T98 


NGC4254 


9.9 ±0.3 


Y95 


1.3 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC4321 


9.2 ±0.1 


H03 


0.4 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC4450 


8.9 ±0.3 


Y95 


-0.2 ±0.2 


CIO 


8.9 ±0.4 


CIO; MAIO 


NGC4536 


8.6 ±0.3 


Y95 


0.3 ±0.2 


CIO 


8.6 ±0.4 


MOlO 


NGC4569 


8.8 ± 0.1 


H03 


-0.1 ±0.2 


CIO 


8.9 ±0.4 


E08; CIO 


NGC4579 


9.0 ±0.3 


Y95 


0.2 ±0.2 


CIO 


9.0 ±0.4 


CIO; MAIO 


NGC4631 


8.5 ± 0.3 


Y95 


0.4 ±0.2 


CIO 


8.4 ±0.3 


MOlO 


NGC4725 


9.1 ± 0.3 


Y95 


0.0 ±0.2 


CIO 


8.7 ±0.4 


MOlO 


NGC4736 


7.9 ± 0.1 


L09 


-0.4 ±0.2 


CIO 


8.7 ±0.4 


MOlO 


NGC4826 


8.1 ± 0.1 


H03 


-0.5 ±0.2 


CIO 


8.9 ±0.3 


MOlO 


NGC5033 


9.3 ± 0.1 


H03 


0.1 ±0.3 


K03 


8.7 ±0.3 


MOlO 


NGC5055 


8.9 ±0.1 


L09 


0.1 ±0.2 


CIO 


8.9 ±0.4 


MOlO 


NGC5194 


9.2 ± 0.1 


H03 


0.4 ±0.2 


CIO 


9.0 ±0.3 


MOlO 


NGC5236 


8.9 ± 0.1 


Y95 


0.4 ±0.2 


CIO 


9.0 ±0.3 


MOlO 


NGC5713 


9.1 ± 0.3 


Y95 


0.6 ±0.2 


CIO 


8.7 ±0.4 


MOlO 


NGC5866 


8.1 ± 0.3 


Y95 


-0.6 ±0.2 


CIO 


8.7 ±0.4 


CIO; MAIO 


NGC5953 


9.0 ±0.3 


Y95 


0.4 ±0.2 


CIO 


8.7 ±0.4 


E08; CIO 


NGC6946 


8.8 ± 0.1 


L09 


0.5 ±0.2 


CIO 


8.8 ±0.3 


MOlO 


NGC7331 


9.0 ±0.1 


L09 


0.4 ±0.2 


CIO 


8.8 ±0.3 


MOlO 



Note. — CIO = ICalze tti ct al.[ (12010); D84 = 'Dufour! (|1984!); EOS = lEnge lbracht ct al.l f2003): F08 = IPukui ct al." J2008J); H03 = 
IHelfer et all II2003I) ; H04 = Hcycr ct al. (2004); H09 = Harris & Zaritsky (2009); 197 = Israel (1997); K03 = Konnicutt ct al. (2003); L79 

= ILequeux et al. l (1 19791 '): L03 = Lcc ct al. (2003 ); L06 = .Lcrov ct al., (,200&'): L09 = iLcrov ct al.. C2009); M97 = .Martin. 11997.); M06 = 
[Mizuno ct al. ('2006); MAIO =IMarble et al.N2010|); MO l O = | Moustakas et al.l 120101) R08 = IRosolowsky fc SimonI ( 120081) : T98 = ITaylor et all 

(1998); W04 = .Wilke et al.l (120041 ): Y95 = IYoung eTal] ( 1T995I) 



alent to 



Mco 



"CO 



^co 



K km s ^ pc^ 



(19) 



with aco = 4.4. Thus, we effectively assume a Milky 
Way conversion factor for CO-emitting gas. It is impor- 
tant to note that our X factor represents the conversion 
from CO luminosity to mass of gas where the carbon is 
predominantly CO, and not the conversion from CO lu- 
minosity to total mass of gas where the hydrogen is H2. 
These concepts are often not clearly distinguished in the 



literature. The approximately constant conversion factor 
derived from virial mass measurements of ex tragalactic 
clouds (jBlitz et al.ll20Q7l: iBolatto et al.ll2008l) (which in- 
clude some but not all of the H2 that is not associated 
with CO) motivates this assumption, though there may 
still be changes in aco of CO-emitting at the factor of 2 
level across the range of metallicities that we study. The 
sense of these would be to increase aco , and to decrease 
the SFR-to-Mco ratio, moving points down in Figure [SJ 
Regardless of systematic effects, the y-axis in Figure [5] 
should be very close to the ratio of ionizing photon rate 
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0.1 1.0 

z 

Fig. 5. — Predicted gaseous specific star formation rate 

SSFR(co,H2, total) = S*/S(CO,H2, total) = ■'^t /*^(CO ,H2 .total) aS 

a function metallicity. As in Figure U the red lines indicate the 
specific star formation rate considering the total gas mass, green 
lines indicate the specific star formation rate for H2, while blue 
lines indicate the specific star formation rate for gas where the car- 
bon is CO rather than C 11. Dashed lines correspond to a total gas 
surface density Stotal = 5 Mq pc~^, dotted lines to Stotal = 20 
Mq pc~^, and solid lines to Stotal = 80 Mq pc~^ . Note that 
these are all total gas surface densities, rather than surface den- 
sities of H2 or CO as on the x axis in Figure |4] Also note that 
we only show values for SSFR^qq) s-i^d SSFR(H2) metallicities 
and surface densities that are high enough for/jjj and /cq to be 
non-zero. 

to CO luminosity. 

To convert measured oxygen abun dances to metallici- 
ties relative to Milky Way we adopt HCaffau et al.ll2008D 

log Z' = [12 + log(0/H)] - 8.76. (20) 

As emphasized by iKewlev fc EllisonI (|200 8') and 
iMoustakas et al.l ()2010D . the adopted calibration has a 
large influence on the metallicities derived from mea- 
surements of strong optical lines. Therefore the overall 
normalization of the metallicities in Figure [5] must be 
considered uncertain by at least ^ 0.1 — 0.2 dex, though 
the internal ordering is likely to be relatively robust. 

A final complication is that most of the literature data 
we have gathered consists of galaxy-integrated values, 
rather than local values. We therefore do not have sur- 
face densities, and we must adopt characteristic values 
in order to compute the relationship between total gas, 
molecular gas, and CO luminosity. Fortunately, the re- 
sults for the specific star formation rate for CO are not 
particularly sensitive to this assumption, as shown below. 

We plot the literature data against our predictions in 
Figure El As the plot shows, there is a clear a correla- 
tion between SSFRco and metallicity that agrees within 
the uncertainties with what one expects if star formation 
follows H2 rather than total gas or CO. Though the data 
are sparse with significant scatter, the difference appears 
to be more than 2 orders of magnitude in the SFR-to- 
CO ratio over about an order of magnitude in metallicity. 
The shape of this trend agrees well with our theoretical 
predictions. Given the uncertainties in estimating several 
of the parameters shown here, we believe this constitutes 
a reasonable first-check on the model. Future improve- 
ments to the measured dust-to-gas ratios, improved esti- 
mates of aco from CO emitting gas, and observations of 



CO from larger samples of low-metallicity galaxies will 
improve the accuracy of the observed trend and allow 
more stringent tests. 

The elevated value of SSFRco in dwarf galaxies and 
the outer parts of spirals relative to the inner parts of 
spirals has been pointed out before (lYoung et al.llTQQGl : 
ILerov et al.ll2006L l2007bHGardan et al.ll2007D . but it was 
unclear if this was due to change in the star forma- 
tion process or a change in the CO to H2 ratio com- 
bined with star formation following H2 rather than CO. 
[Pelupcssy & Papadopoulos (2009) argued that the ele- 
vated value of SSFRco could be explained if the latter 
were true, but they simply assumed that star formation 
was correlated with H2, rather than explaining the cor- 
relation from first principles. Here we have provided the 
missing explanation, and Figure [5] demonstrates that a 
model based on it can quantitatively reproduce the ob- 
servations. 

4.3. Predictions for Future Observations 

Figures S] and [5] also contain clear predictions for ob- 
servations. At present it is extremely difficult to measure 
either total gas masses or H2 masses in galaxies with low 
metallicity, because CO breaks down as a tracer of H2 in 
these environments, and because the H2 itself does not 
emit. Observations are available for only a few galaxies 
based on using pro xies other than CO for the H2 (e.g. 
iLerov et"al]|2007al) . However, future dust observations 
with facilities such as Herschel and ALMA will make it 
easier to obtain H2 masses, and thus total gas masses, 
for more galaxies. Our work contains a clear prediction: 
for these galaxies, the specific star formation rate with 
respect to H2 mass should be essentially independent of 
metallicity, while the specific star formation rate for the 
total gas mass will behave in the opposite sense as for 
CO: low metallicity galaxies will have lower SSFRtotai, 
even as they have higher SSFRco- 

With resolved observations an additional test becomes 
possible. Figure m shows that SSFRtotai is essentially flat 
at high surface densities, but turns down sharply at sur- 
face densities below a metallicity-dependent value. Such 
a drop is seen below Stotai ^ 10 M(d pc~^ in observations 
of Solar metallicity galaxies IBigiel et al.l (2008). Mod- 
els in which the star formation rate depends on global 
gravitational instability or similar phenomena that do 
not care about gas cooling or chemistry predict that 
the surface density at which the total gas specific star 
form ation rate dro ps should not vary with metallicity 
(e.g. lLi et al]|20060 . In contrast, our work here suggests 
that it should scale roughly inversely with metallicity; 
ISchavd (l200l . using a thermally-based model that an- 
ticipates some of this work, makes a similar prediction 
(his equation 25), although the metallicity-dependence 
in his model is weaker than in ours. Resolved measure- 
ments of the total gas surface density in nearby galaxies 
should be able to test which prediction is correct: does 
SSFRtotai always change sharply at ~ 10 Mq pc~^, or 
is that value metallicity-dependent? Preliminary results 
indicate appear consistent with a dependence on metal- 
licity (Fumaga lh et al...2010.) , but a more systematic sur- 
vey is needed. 

5. SUMMARY 
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In this paper we consider which phase of the interstel- 
lar medium should correlate best with the star formation 
rate, and why. Our wor k extend s the p revious examina- 
tion of this problem by iSchavd ()2004[ ) . We use chem- 
ical and thermal models of interstellar clouds to show 
that the star formation rate is expected to correlate most 
closely with the molecular hydrogen content of a galaxy. 
In order to undergo runaway gravitational collapse to 
form stars, gas must be able to reach low temperatures 
and therefore low Bonnor-Ebert masses. Its ability to do 
so, however, is impaired by the interstellar radiation field, 
which heats the gas. Only in regions where the ISRF is 
sufficiently excluded by extinction can star formation oc- 
cur. We find that such regions are also the regions where 
the gas is expected to be predominantly H2 rather than 
H I, and for this reason the star formation rate correlates 
with, but is not caused by, the H i to H2 transition. In 
contrast, the chemical makeup has relatively little effect 
on the ability of the gas to cool. All of these results are 
robust against a very wide range of variation in metal- 
licity, radiation field, or other properties of the galactic 
environment. 

That star formation correlates with II2 rather than ei- 
ther total gas mass or CO mass has strong obseryationa l 
implications (see also iPelupessv fc PapadopoulosI 120091 ). 
The fraction of H2 gas where the carbon is in the form 
of C II rather than CO is a strong function of metallicity. 
Galaxies with low metallicity tend to have large masses 



of H2 where there is fairly little CO. If this material is 
able to form stars, as we predict, then the star formation 
rate per unit CO mass should be very large in low metal- 
licity galaxies. We see that precisely this phenomenon 
is found in observed galaxies. Finally, we note that the 
fraction of the total gas mass where the hydrogen is II2 
rather than H i is also a strongly increasing function of 
metallicity. Since star formation correlates with II2, we 
predict that the star formation rate per unit total gas 
mass should be small in low metallicity galaxies. This 
prediction can be used to test our calculations. 
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APPENDIX 
A. PARAMETER CHOICES 

All molecular and atomic information, including level e nergies, Einstein coe fficients, and collision rate coefficients, 
is taken from the Leiden Atomic and Molecular Database fSc hoier et al.|[2005l ). Other input parameters are gathered 
from a variety of sources, and our fiducial values and references are given in Table[TJ Most of these are straightforward, 
and where a quantity has been observed in the Milky Way, we extrapolate to other galaxies by assuming that element 
abundances and dust to gas ratios are simply proportional to metallicity. Parameters that are not directly observed 
or extrapolated we discuss in the remainder of this Appendix. 

A.l. Dust-Gas Thermal Exchange Rate Coefficients 

For the dust-gas thermal exchange rate coeflncient ad sr, we adopt a standard value of 3.8 x 10"'^^ erg cm^ K~^/^ in 
Il2-dominated gas in the Milky Way (IGoldsmithll200l . and we assume that this will also hold in galaxies of similar 
metallicity. We must extrapolate this both to lower metallicity galaxies, and to predominantly H i gas. For the former, 
we assume that the total surface area of grains is proportional to the metal abundance. For the latter, we must 
account for both the change in the number of particles and the masses of the individual particles, which alters their 
speed. We assume that the accommodation coeffi cient is equal for H, Hg, and He. (More accurate approximations are 
possible if the grain size distribution is known, e.g. lHoUenbach fc McKed ()1979f ). but given the uncertainties in how size 
distributions vary from galaxy to galaxy and the relative lack of importance of this effect, we omit this complication.) 
With this assumption, the dust-gas energy exchange rate provided by gaseous members of species X is proportional 
1/2 

to nx/tn^ , where ux is the number density of members of that species and mx is the mass of that species. Thus in 
regions of atomic and molecular hydrogen respectively we have 

"HI , "Ho / . , ^ 

ade,moc—^ -\ -7^ (Al) 

"H2 , "He 
"^H2 "^Hc 

where the constants of proportionality are the same in each case. In a predominantly atomic region njji — "h, while in 
a predominantly molecular region = riH/2; in both cases rifjo ~ "h/10. Similarly, mfjc = 4mH and = 2toh- 
Combining the dependence on metallicity with that on chemical phase, we arrive at our final expressions for adg- 

adg32 =3.8 X IQ-^^Z' erg cm^ K"^/^ (A3) 
adg,Hi = 1-0 X lO^^^Z' erg cm'"^ K'^^^. (A4) 
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For the radiation temperature, Tj-ad, as discussed above we must select a single value to characterize the re-radiated 
infrared field within a cloud. We adopt Trad = 8 K for this, near the minimum temperature seen in ammonia 
observations of Galactic cold clouds (iJiiina et al.lll999l) . Both this choice and our choice of adg have very little impact 
on our results due to the weak dust-gas coupling in the density range with which we are concerned. 

A. 2. Opacities 

Our calculations depend on the absorption cross section per H nucleus aa to UV photons, and the visual extinction 
per unit hydrogen column Ay /N-a. Note that aj, reflects only absorption, while Ay includes scattering as well. These 
quantities depend on the extinction curve and vary between dense and diffuse environments even within a single 
galaxy at a single metallicity, and thus t heir exact values are uncertain by a factor of two. As a guide we examine the 
values given by the models of Wcing artner fc Drainel (j2001f) for the Milky Way, Large Magellanic Cloud, and Small 
Magellanic Cloud0 Once we scale the LMC and SMC models to the same total dust to gas ratio as the Milky Way 
models, we obtain valu es of (7^/10"^^ cm'^ = 0.7, 1.7, 2.0, 1.8 and {Ay / Nn) / 10~'^^ mag cm'^ = 4.5, 5.1, 3.6, 2.8 for 
iWeingartner fc Drainef s models for Milky Way sightlines with Ry = 5.5, Milky Way sightlines with Ry = 3.1, the 
average LMC sightline, and a sightline through the SMC bar, respectively. We therefore adopt intermediate values 
ad = 1-0 X cm^ and Ay/N^ ~ 4.0 x lO^^^Z' mag cm^. With these fiducial choices, the UV absorption is 

larger than the the visual extinction by a factor of 2.5/1.08 = 2.31 (where the factor of 1.08 accounts for the conversion 
from magnitudes to true dimensionless units). 

A. 3. Far Ultraviolet Radiation Field 

The FUV radiation field, parameterized by Gq, affects both gas temperature and chemistry. There is no unique 
value of Gq that characterizes all gas clouds across all galaxy types, and even within a single star-forming cloud the 
ambient UV field increases with time as stars form. As a fiducial choice we adopt Gq = 1, the Solar neighborhood 
value. In regions where star formation is ongoing the mean is likely larger, but we are interested in the initiation of 
star formation, so it seems prudent to select a value typical of regions before the onset of star formation. 

To test our sensitivity to this choice, we have recomputed our model grids with values of Gq — 0.1 and Gq ~ 10, 
and re-plotted Figures [1] - [3] of the main text in Figure [T] Not surprisingly, the H i to H2, C ii to CO transitions, 
and warm to cold gas transitions all move to higher density and Ay for higher Gq, and to lower density and Ay for 
lower Gq. Nonetheless, over a factor of 100 range in Gq, we will retain the excellent correlation between the H i to 
H2 transition and the drop in temperature from hundreds of K to ~ 10 K, along with the concomitant drop in the 
Bonnor-Ebert mass from several 1000 Mq or more to a few Mq. At all three values of Gq, contours of constant 
align closely with contours of constant Bonnor-Ebert mass. In the lower two panels, note that, contours of constant 
Mbe remain largely vertical, indicating that the Bonnor-Ebert mass drops systematically as the II2 fraction increases. 
Thus our conclusions are robust against large variations in Gq. 

A. 4. Cosmic Rays 

The cosmic ray heating rate is uncert ain in two ways. First , the energ y yield per prim ary cosmic ray ionization qcn 
is 6.5 eV in predominantly neutral H i ([Dalgarno fc McCravl[r972: Wolfire et al.l 119951) . In H2 the yield is uncertain. 
Both dissociative recombination of H2 and excitation of its rotational and vibrational levels followed by collisional 
de-excitation provide extra channels for energy transfer from primary cosmic ray electrons to thermal motion, but the 
importance of these processes is likely density-dependent. Values given in the literature range frorn near ly the same 
energy yield as in H i up to 20 eV per primary ionization (iGlassgold fc Langeilll973l: iDalgarno et al.|[l999l ). We follow 
I Wolfire et al.l ()2010f ) in adopting an intermediate value qcR = 12.25 eV, but this should be regarded as uncertain by a 
factor of 2. 

A similar uncerta inty affects the cosmic ray ionization rate. For Milky Way-like galaxies we adopt C = 2 x 10^^^ s~^ 
(| Wolfire et al.llMfll ). but the p rimary cosmic ray io nization rate even in the Milky Way is substantially uncertain, and 
may be hig her than this value (|Neufeld et al.l[20To[ ). The cosmic ray intensity also varies with the star formation rate 
in galaxies (lAbdo et al.ll20ldal) . Thus in low metallicity galaxies, which also tend to have low star formation rates, the 
cosmic ray ionization rate is likely lower than the Milky Way value. The scaling is extremely uncertain; for lack of a 
better alternative, we simply take the cosmic ray ionization rate to scale with the metallicity. 

To test our sensitivity to our choice of cosmic ray ionization rate, we have recomputed our temperature grid with 
cosmic ray ionization rates that are a factor of 3 larger and a factor of 3 smaller than our fiducial choice, thereby 
exploring a ~ 1 decade range in cosmic ray heating rate. We plot the results in Figure [21 Again, we see that the 
results are robust, in the sense that contours of constant Bonnor-Ebert mass remain closely aligned with contours of 
constant H2 fraction as we vary the cosmic ray ionization rate. 

We do caution, however, that this breaks down if we select a cosmic ray ionization rate that is extremely high (more 
than 10 — 30 times our fiducial one). For such high cosmic ray ionization rates, cosmic ray heating becomes more 
important than grain photoelectric heating even in unshielded, low Ay gas. In this case the temperature and Bonnor- 
Ebert mass become uncorrelated with Ay, and vary with density alone. However, such high cosmic ray ionization 
rates are inconsistent with observations in the Milky Way. 

* IDraind 120031) gives updated models for the Milky Way, but 
since these are not available for the LMC and SMC, we use the 
older models instead. The difference is at most tens of percent. 
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Fig. 1. — Left column: same as Figures [T1 l2l and|3] but with Gg = 0.1 instead of Gq = 1. Right column: same as Figures [T] [J] andlSl 



but with Gq = 10 instead of Gq = 1. 
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